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We present a study of the interaction of a relativistically planar channeled electron with an intense 
electromagnetic field. Using a S-Matrix approach in the Strong Field Approximation, it is shown 
that the crystal periodicity affects drastically the excitation process, suppressing the possibility 
of multiphoton absorption except for some particular cases. This selective excitation opens the 
possibility to control the dynamics of the channeling process by means of an external field. Explicit 
expressions for the S-matrix N-photon excitation rates together with the corresponding conservation 
laws are obtained from the relativistic quantum mechanical Dirac equation. 



PACS: 61.85.+p, 03.65. Pm, 12.20.Ds 



I. INTRODUCTION 

Channeling in crystal lattices occurs when an acceler- 
ated charged particle is introduced into a crystalline tar- 
get at sufficiently large energy. Depending on the crystal 
orientation, the particle's trajectory may be aligned with 
a major crystal direction and the penetration may reach 
anomalous depths. Although the possibility of this ef- 
fect was already pointed out very early by Stark |Q, it 
was demonstrated experimentally 50 years later by Rol 
et al. (2| , when the result of the ion sputtering was found 
to depend strongly on the orientation of the target crys- 
tal. After the discovery, the theoretical and experimental 
work increased rapidly and extended to the case of chan- 
neling of electrons and positrons (3) . 

In this paper we will investigate the excitation dynam- 
ics of a planar channeled electron under the influence 
of an external electromagnetic field. As a main result, 
we demonstrate that the crystal periodicity introduces a 
momentum conservation condition which affects the ef- 
ficiency of the different N-photon channels of excitation, 
leading to the strong suppression of photon absorption 
in a broad range of situations. The multiphoton exci- 
tation of channeled particles has been already addressed 
by Avetissian et al. iQjHJ by assuming an electromagnetic 
wave copropagating with the electron, and with a fre- 
quency resonant to the (Doppler-shifted) lower-energy 
level transitions. In the present case, however, we are 
interested in a complementary situation where the chan- 
neled electron is excited to a final state lying in the crys- 
tal quasi-continuum. Since the transition is produced by 
the interaction of the electron with an external intense 
optical field, the strong field approximation (SFA) con- 
stitutes a more appropriated procedure in comparison to 
the discrete level approach in j| . 

SFA theories have been developed in the context of 
ionization of atoms in strong laser fields. Among them, 
the so-called Keldish-Faisal-Reiss (KFR) theory [0-0 is 



based on the S-matrix approach, where the final state 
is approximated by a Volkov state, which describes the 
evolution of a free electron driven by the electromag- 
netic wave. SFA theories describe most of the relevant 
aspects of the atomic ionization, including multiphoton 
absorption and multiphoton excitation above the ioniza- 
tion threshold (ATI). 

Although employed mainly in the atomic and molecu- 
lar context, S-matrix SFA approaches can be used in any 
general situation in which the field interaction energy is 
comparable with the energies of the matter system. In 
fact, for the higher energy boundstates, the intensity of 
the field required to promote an electron to the contin- 
uum does not have to be very high, and yet SFA can be 
used. On the other hand, SFA requires the matter po- 
tential to be approximately constant over the complete 
interaction time. In our case it suffices with a moderate 
intensity field, (10 12 — 10 13 W/cm 2 ), while the crystal sta- 
bility can be ensured by a sufficiently short pulse (about 
100/s), which still enclose enough cycles to ensure the 
adiabatic limit involved in the theoretical approach. It 
should be mention that, in the case of atom ionization by 
strong field, the adiabatic assumption is correct even in 
the case of a few cycles pulse-length, and there is no rea- 
son to think that in this channeling case the thing should 
be different. 



II. GEOMETRY OF THE SYSTEM AND 
DESCRIPTION OF THE UNPERTURBED 
CHANNELED ELECTRON. 



Let us consider the interacting geometry depicted in 
fig. [I]. A relativistic electron, with velocity parallel to 
the x — axis, is introduced in a crystal while interacting 
with a copropagating or a counterpropagating electro- 
magnetic wave. For particular crystal orientations, the 
electron is confined transversely to trajectories close to 
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the initial injection axis. Our case of planar channeling 
occurs when the electron is injected parallel to a crystal 
plane |||| . We will assume an electromagnetic planewave 
field, which propagates in the same direction of the in- 
jected electron, and linearly polarized in the y — axis, i.e. 
orthogonal to the crystal plane. 

Since the channeled particle is injected with a relativis- 
tic velocity, the crystal potential may be well approxi- 
mated by a spatial average over the crystal plane coor- 
dinates, as it is done in the so-called continuum model 



V(y) 



W 2 r L z /2 



L X L Z J ~L m /2 J -L z /2 



U(x,y, z)dxdz (1) 



where L x and L y are the crystal plane dimensions, and 
U(x, y, z) is the crystal potential. With this effective po- 
tential the electron motion in the crystal channel, driven 
by the external e.m. field, will be confined to the polar- 
ization plane xy. 

In the most general case, the quantum description of 
the electron's dynamics is described by the Dirac equa- 
tion 



{ca • (p - eA/c) + /3mc 2 + V(y)} $(x) = £ B $(x) (2) 

Let us first consider the unperturbed channeling situ- 
ation, A = 0. Since the averaged potential depends only 
on the y coordinate, a general positive energy solution of 
the channeled electron can be written as 



*(X) = / d P y 



mc 
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Where u px Py is the positive energy solution of a free 
Dirac electron 10 llj], E Pj . :Py = ^Jc 2 p 2 + c 2 p 2 + m 2 c 4 , 

/ 1 \ 



Px,Py 



2mc 2 
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(4) 



Introducing (||) into (||), with A = 0, one can obtain 
the following form for the Dirac equation in momentum 
space: 



E BUp xiP € Pv = {ca xPx + ca yPy + /3mc 2 }«J tiP> ^ + f dp' y J P °" Py V Py - p > y u PxiP ' y £ p , y (5) 



Being V p the Fourier transform of the interplanar poten- 
tial at the spatial frequency p/h. Due to the nature of the 
averaged potential, the channeling along a crystal plane 
is only possible if the energy of the electron's transversal 
dynamics is moderate, i.e. non relativistic. In this case, 
it is justified to approximate equation (|J) to second order 
of cpyjE Px . 
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By introducing eq. (||-||) into (jE)), and since the scalar 
potential V Py - Py is a first order term in cp y /E x , the Dirac 
equation is reduced to the same identity for every non- 
zero component of the spinor u px 



>Pv' 



E B t Pv = + E Px }£ Py + I dp' v V Pv - p , y ^ ly (9) 

Px J 

Computing the inverse Fourier transform in the y coor- 
dinate, we finally end with a Schrodinger-like equation 
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Note that, for the field intensities considered here, this 
approximation will remain equally valid when considering 
the interaction with the external electromagnetic wave. 



2j x m 



v(y)H(y) 



(10) 



where E P7 . — ^ x mc 2 , and ts = Eb — j x mc 2 corresponds 
to the non-relativistic eigenstate energy. 

Several explicit forms of the averaged crystal potential 
may be found in the literature ||. Among them, those 
derived from the Thomas-Fermi screened two-body po- 
tentials have been widely used j^^], and have been fine 
adjusted by standard Hartree-Fock many-body calcula- 
tions [|l4j]. From the theoretical point of view, model 
potentials are more convenient, since they allow for fur- 
ther analytical work while keeping the essential features 
of the interaction. For instance, V(y) = —4Voy 2 /d 2 was 
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proposed by Avetissian et al. jj| in the context of com- 
putation of the multiphoton transitions between bound 
states of the continuum potential, for a positron inter- 
acting with a strong electromagnetic wave. The form 
V(y) — — Vq cosh~ 2 (y/6) has also been used to study the 
possible amplification of x-ray channeling radiation 
This later form has a better resemblance with the av- 
eraged Thomas-Fermi potentials while still allowing for 
an analytical diagonalization and, therefore, we shall use 
this potential for our calculations. The transverse energy 
spectrum for this case can be cast in the following form 



2b 2 m"/ s 



■(s-nf 



(11) 



where 



can be 0, l,...,[s], being s = 



1 + 2& 2 m7 x y _ Note ttxatj ag a regult of tne spat j a i 

averaging, the scattering with the continuum potential 
affects only to the transverse dynamics of the channeled 
particle, while the very large longitudinal momentum re- 
mains unaffected. 



III. S MATRIX DESCRIPTION OF A 
CHANNELED PARTICLE INTERACTING WITH 
AN ELECTROMAGNETIC WAVE. 



can be exported to any other system, provided its eigen- 
states can be found analytically. To our knowledge, how- 
ever, this is the first time that it is applied to the rel- 
ativistic channeled electron in interaction with an elec- 
tromagnetic wave. In the present case, $j corresponds 
to the unperturbed channeled electron state discussed in 
sec. ||, and ^P* is an arbitrary final state, solution of 
the complete equation (||) . Since this exact solution is not 
available, the success of the S-matrix approach consist in 
finding a suitable approximation. In the strong field ap- 
proximation, the interaction with electromagnetic field 
is assumed to be the relevant for the final state, there- 
fore ^ is approximated in terms of the Volkov states, 
\&V"(x» t) |p2] , p3| . These wavefunctions are solutions of 
eq. (||) for V(y) = and A(x, t) ^ 0, and describe a free 
electron in the presence of an electromagnetic field. The 
form of these states for a laser field pulse of arbitrary 
form is 1231 : 
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Let us now add the electromagnetic excitation to the 
problem. As it is well known, eq. (|^) does not accept 
analytical solutions for a space-time dependent vector 
potential. In such situations, the S-matrix approach of- 
fers a standard procedure to find approximated solutions 
fl5f , used specially in quantum field theory (l(] |l(],[ljj and 
scattering |l8|. 



The S-matrix approach takes as a starting point the 
following exact relation 

* f (x)=9 v (x) + J d 4 i'Gy(i,i')7^(s')*/(^') (15) 

where Gv(x,x') is the Volkov Green's function 0], and 
keeps the lowest order term in powers of V(x). The tran- 
sition amplitude obtained is pO|,pl[: 



A. The general relativistic case 

The relativistic SFA S-matrix theory for the Dirac elec- 
tron in an atom can be found in |Ts|-|2l|] . The general ex- 
pression for the transition amplitude using time-reversed 
S matrix theory has the following form 



Sn = lim 



(12) 



Although mainly used in the strong field ionization of 
atoms and molecules, this approach is quite general and 



hc \ (2tt) 3 E 



d^xe~ iS vL 



(eA ^ 2 {A-A){ \ 
\ e4 ' 2c(s.p) ) 



(16) 



where we have used ( |l3| ) and we have assumed the trans- 
verse character of the e.m. field, s ■ A = — s ■ A = 0. 
Separating the time and the space integrals in equation 
(Ea), the transition amplitude takes the form: 
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being and ^S2(t): 
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Once the transition amplitude (S — l)fjj is known, 
the total transition rate can be computed as: 



W = 



Vdp 

(2^)3 



(20) 



where is the normalization volume and w is the tran- 
sition probability per unit of time, which is defined as 



1 



w 



lim -\(S-ir f 



(21) 



Another important magnitude is the transition rate per 
unit of solid angle, which has the form: 



dW 
~d£l 



V 



(2tt) s 



wp 2 dp 



(22) 



It should be mention that SFA S-matrix theory cannot 
be considered as a perturbation series in powers of V(x) 



since the initial state <f>i is an eigenstate of the potential 
itself including all the properties of the crystal. The pres- 
ence of this wavefunction produces a new behavior in the 
scattering section, introducing all the differences with the 
atom ionization or the simple Compton scattering. 



B. Application to the case of a monochromatic and 
linearly polarized laser field 

Lets us focus our attention to the geometry depicted in 
fig. 0, where the electromagnetic field can be a copropa- 
gating or counterpropagating linear polarized planewave 
of frequency u, A.(<p) = Aj(tp)ej = Aq cos(wi — k ■ x)ej. 
where e,- is the field's polarization vector. 



The phase factor of the Volkov function ( |l4|) now reads 



as 



r/ V e 2 A\U l \ epjAo . .. . r e 2 A\ . , nn 

S = 6 \ k~ + AheWTpj ) *» + h^tpj Sm(fc ' X) ~ e Sm(2(fc ' X)) (23) 

The resulting exponential factor can be expanded as a series of Bessel functions 

e -iS =e ^ r [^+j^)\^ JN+2n(r])Jn(0e- iN{k - x) (24) 

N,n— — oo 

where rj and £ are factors which only depend on the momentum of the Volkov function, and the frequency and 
amplitude of the laser field: 

= epjAp = e r e 2 A 2 

1 hc(k-p) ? &hc 2 (k-p) v ' 

Substituting (H) in eqs. (jl|) and @: 

= - L ir E JN +2n (v)Jn(0e* — s — Mm) 26 

h e - — ' 77 

N,n=-oo 



and 



= ~ j dxe- lS |A| 2 e- 1 -^^ 

= — Z- 2^ 1 2 JN+2n(V) + ^-(^JV+2n+l(?7) ~ JjV+2n-lWJ 

xJ„(C)e' $ a ( qjv ) (27) 
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Where q^, u>n and 5>,(qjv) are denned as: 

UN 



The time integrals appearing in (|17| ) can now be calculated as 
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(31) 



and 



/ = E 



N.n= — oo 
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being T = t — to 

To compute the rate of excitation, we should use this two equations together with eq. ( |T7j ) to calculate 



with 
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(35) 
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Substituting the time integrals, [31] and |3^, in each term, we may calculate the transition probability per unit of 
time as: 



1 



w= lim -|(S-1)^T= lim 



where 
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Finally, the excitation rate is 



in 2 E\ h J 1 1 c{k-p) * ' y " \2c(k-p) 



2 



^) £ - ^ 2 Re(A 2 ) + f-^ ) StAA 5 ^ N - 2fe) (43) 



Where we have defined: 



Sl= £ (JV + 2n) J j V+2nWJn(0 (44) 
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»? 
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n— — oo V ^7 ^ / 

Ai = $+( qj v) 7 J+ 7 ">p + 7V^(qjv) (47) 
A 2 = $+(qAr) 7 - ? ' + 7 ^>p + 7 ^(qiv) (48) 
A 4 = $+(q JV )^ + 7 u pWp + 7 ^(q J v) (49) 



IV. CONSERVATION LAWS AND CLOSING OF 
EXCITATION CHANNELS. 

As expressed in eq. (|43|), the transition probability 
is a function of the initial momentum-space probability 
amplitude, $j(qj\r). Let us now assume injected electron 
of positive energy with a wavefunction of the form (|3|), 
therefore e r = +1. From the delta function in eq. ( |43[ ), 
we obtain the following energy conservation relation: 



4hc 2 (k -p) 



N 



hw = E 



B 



(50) 



On the other hand, an additional conservation law relates 
the momentum of the final and initial states, p and p; 
respectively, in (p8|). For a positive energy electron, this 
reads as 



Pi = T 



Ahc 2 {k -p) 



N 



(51) 



Since the electromagnetic field propagates along the x- 
axis, this condition may be splitted into two parts 



(Px)i 

(Py)i 



Px ± 
Py 



Ahc 2 (k -p) 



- N 



hk 



(52) 
(53) 



with (k ■ p) — fc(-^ TPx): an< 4 where k = |k|, the top 
sign refers to a field copropagating with the electron, and 
the bottom to the counterpropagating case. These three 
equations, (|50|), ( |i>2| ) and (|53"|), describe the energy and 



momentum changes due to the stimulated absorption or 
emission of N photons. Combining these with the energy 

expression for the final state, E = yjc 2 p 2 + c 2 p 2 + m 2 c 4 , 
we obtain a closed formula for the energy conservation of 
the multiphoton process, in terms of the initial momen- 
tum and the field parameters 



'■(Px)l +C 2 (Py)1 + 



El 



2"fBlTLC 2 
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where wc have defined the initial energy of the electron as 
Eb = Jbitt-c 2 , and the initial relativistic velocity factor, 



(Px)i = ( Px)i/lBmc. S ince E B 
(7*)< = VV(1 - (Px)'f) we have 



es + (lx)imc 2 , with 



M 
2m 



A 2 e 2 

(7«)iefl + (55) 
4mc 



The interpretation of this energy conservation relation is 
straight-forward if we take as reference system a frame 
propagating with the electron, with its initial velocity 
(Px)i/(jx) l m. The frequency u>' = w 7s (1 =F cor- 
responds to the Doppler-shifted electromagnetic wave, 
and (j x )i£B is the result of the Lorentz transform of the 
bound state energy ( pd|) from the laboratory to the mov- 
ing frame. Equation (|55|) states the resonance condition 
for a (Doppler-shifted) TV-photon transition from a bound 
state to a state lying in the crystal quasi-continuum of 
momentum given by eqs. (p2) and (J53|). This is, in 
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essence, the crystal equivalent to the ionization of atoms 
by intense fields. Note, however, that in the atom case 
a ionization channel for any photon number N is always 
possible, since the initial state is distributed continuously 
over the momentum space and, therefore, a non-zero 
transition probability exists for any (p y )i which fulfills 
the condition (|55|). This is not the case for the chan- 
neled electron, since the crystal plane periodicity forces 
a discretization of the electron states in the transverse 



energy of the electron: 



coordinate of the momentum space p y 



~ 2-Kh 



nAp y , 



being n an integer and d p the interplanar distance. As a 
consequence, in the general case the TV-photon channel of 
excitation should be strongly suppressed, except in those 
particular cases in which 



N ■ 



{lx)j£B 
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n 2 Ap y 
2m 
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Xmc 2 ) 



(56) 



holds for TV and h as integer numbers. As a consequence, 
this opens the possibility of selective excitation of chan- 
neled electrons in terms of their initial velocity, or per- 
mits its control through the variation of the electromag- 
netic field parameters. Figure ^a-d show the possible 
A-photon channel excitations as a function of the initial 
electron's energy and for the lowest orders of transverse 
momentum transferred n. Each plot shows the result 
for a different initial channeling bound-state. We assume 
planar channeling along the (110) plane of Si by selecting 
the potential parameters Vq = 20.4 eV and b = 0.03 nm 
reproducing [ ^5| , and a counterpropagating TiSa laser of 
3.51 x 10 12 W/cm 2 (A ~ 800nm). Note that the number 
of photons A should be an integer quantity, therefore 
the figure shows clearly that, except for very particu- 
lar choices of the electron's initial energy, the excitation 
channels are closed. 

The same figure can be done for the case of a copropa- 
gating electromagnetic field. Figure ||a-d show again the 
possible A-photon channel excitations as a function of 
the initial electron's energy, for the lowest orders of trans- 
verse momentum transferred assuming the same crystal 
and laser parameters as in figure |^. An important in- 
crease of the number of photons needed to excite the 
electron, attributable to the Doppler redshift, can be ob- 
served. Under this circumstances, even when the en- 
ergy and momentum constrains are fulfilled, the process 
may involve a very small transition probability due to the 
high number of photons needed. To give an idea of the 
order of this probability one can make use of the asymp- 
totic expansion of the Bessel functions for large orders, 

J »W " vfe(^)"' & beingne = lim^oo^ + i)", 
to calculate a limit of the number of photons above which 
the transition probability will be negligible. The criterion 
to be used here will be to consider negligible the Bessel 



=1*1 



< 0.1. Applying it to our case one 



function when 
obtain the limit of the number of photons for each trans- 
verse momentum transferred as a function of the final 
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where 



ary(l — /3 X ). Assuming that the electron 



finishes in a state of the crystal quasi-continuum and 
that the energy along the axial direction do not change 
significantly during the evolution, one can approximate 
7 w (7x)j and (3 X « {(3 x )i. Consequently a/, and there- 
fore NLimit, can be expressed as a function of the initial 
electron's energy. Figure [§a-d show with arrows the point 
when the number of photons required for the excitation 
process surpass the NLimit- For energies above this point, 
the transition probability reduces drastically and we can 
consider that no excitation takes place, even though the 
energy and momentum conservation relations may be ful- 
filled. It should be pointed out that eq. (|57|) is only valid 
for the case of large orders in the Bessel functions. This 
means that it should be taken qualitatively in all cases in 
which NLimit is a small quantity, as for instance in figure 
||a for the case h = 0. Note also that those cases in which 
the arrow is not shown correspond to NLimit outside the 
plotting region, i.e. the A-photon excitation is possible 
along the complete plotted line. 

Finally, let us remark the fact that the photon exci- 
tation number is greater than in the counterpropagating 
case increases the sensitivity of the channel process to 
the selective excitation in terms of the laser parameters. 



V. CONCLUSIONS 

We have computed the explicit forms of the S-matrix 
transition probabilities for the N-photon absorption of a 
relativistic electron channeled along a crystal plane. In 
contrast to previous works, we consider the interaction 
with an intense electromagnetic wave, generated exter- 
nally, which may excite the electron to high-energy states 
lying in the crystal quasi-continuum. Due to the crys- 
tal periodicity, we show that the energy and momentum 
conservation equations constraint strongly this excitation 
process, suppressing the multiphoton absorption except 
for some particular cases. Under these circumstances, the 
selection of a single multiphoton channel of excitation is 
feasible by an adequate choice of the external laser pa- 
rameters, opening a broad range of possibilities for the 
coherent control of the channel electron's dynamics. The 
case of an electromagnetic field copropagating with the 
injected electron is also studied showing an important 
increase of the number of photons needed to excite the 
electron due to the Doppler redshift. For this case, we 
give an estimation of the maximum photon number for 
which the excitation process is not negligible. The se- 
lective excitation in the copropagating case is found to 
be more sensitive to the electron's energy and the trans- 
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verse momentum transferred in the transition than in the 
counterpropagating one. 
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FIG. 1. The system to be studied consists of a relativistic 
electron channeled in a crystal and interacting with a copropa- 
gating (top picture) or a counterpropagating (bottom picture) 
electromagnetic field. 
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FIG. 2. This figures shows the values of (*fx)i needed to 
open a specific N-photon excitation channel for the lowest or- 
ders of transverse momentum transferred h. All the pictures 
correspond to an electron channeled along the (110) plane of 
silicon, driven by a counterpropagating linear polarized TiSa 
laser, A = 800 nm, of 3.51 x 10 12 W/cm 2 . The different pic- 
tures represent distinct initial bound states of the interplanar 
potential. The top one corresponds to the ground state case 
(n = 0), the second picture to the first excited state (n = 1) 
and so on. The continuous, dotted and dashed lines represent 
the excitation process with zero, one and two quanta of trans- 
verse momentum transferred (h — 0, ±1, ±2) respectively. 



FIG. 3. The same situation as in figure |2j but with the ex- 
ternal electromagnetic field copropagating with the channeled 
electron. The arrows represent the Nhimit for each transition 
line. The thinner arrow in figure (a) must be only considered 
as a qualitative estimation (see text). The transition lines 
without arrow mean that the NLimit occurs for parameters 
beyond the plotted region. 
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